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The role of local fields in the optical response of silicon nanocrystals is analyzed using a tight 
binding approach. Our calculations show that, at variance with bulk silicon, local field effects 
dramatically modify the silicon nanocrystal optical response. An explanation is given in terms of 
surface electronic polarization and confirmed by the fair agreement between the tight binding results 
and that of a classical dielectric model. From such a comparison, it emerges that the classical model 
works not only for large but also for very small nanocrystals. Moreover, the dependence on size of 
the optical response is discussed, in particular treating the limit of large size nanocrystals. 
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Silicon nanocrystals (Si- nc) constitute very promising 
building blocks for the next generation of optoelectronic 
devices. Since the discovery of a strong optical activity of 
Si-ncji their fabrication techniques have been much im- 
proved. Nowadays, small (< 2 nm) Si-nc samples with 
narrow size distributions and high quantum yields are 
produced by chemical synthesisi^i&i An open discussion 
in literature is on the influence of the passivant on their 
optical response. Recent experiments have shown that 
new kinds of passivation by organic molecules may open 
new scenarios for enhancing the Si-nc quantum yieldi^ 
The passivation and the background in which Si-nc are 
synthesized cannot therefore simply be ignored, and their 
role still must be carefully examined. From a theoreti- 
cal point of view, some work has been done in analyzing 
the effects of passivant species in silicon nanocrystals. 5 
But very few papers discuss the contribution of the back- 
ground from an atomistic point of view. Indeed, stan- 
dard calculations are based on the independent particle 
approach, and they fully neglect the presence of a back- 
ground embedding medium, treating the nanocrystal as 
an isolated object. In this paper, using a tight bind- 
ing framework, we show that local field effects (LFEs) 6 
deeply influence the Si-nc optical response. We compare 
our results (i) with a dielectric model based on classical 
electrostatics, (ii) with first principles calculations, and 
(iii) with experimental data. In the following, the dielec- 
tric model and the atomistic approach are illustrated. 
The results are then discussed. 

There are two remarkable differences between the op- 
tical response of bulk Si and Si-nc. The first is the 
quantum confinement effect (QCE). Quantum confine- 
ment has given for years a qualitative and often quantita- 
tive description of the Si-nc photoluminescence blueshift 
upon decreasing the nanocrystal size, and it has been 
widely studied using both semiempirical and first princi- 
ples computational toolsi 7 - The second difference, whose 



importance, in our opinion, has always been underesti- 
mated, is the surface polarization effect (SPE). When a 
dielectric structure is embedded into a background char- 
acterized by a different dielectric constant, its optical re- 
sponse can dramatically change. Indeed, as it is well 
known from the classical electrostatics, the charge accu- 
mulation due to the dielectric mismatch across the sur- 
face causes the macroscopic field inside to be very differ- 
ent from the external field£ Numerous definition of the 
dielectric constant of a nanocrystal can be given iZ-uLiii We 
follow a procedure, in which, starting from the polariz- 
ability a of a Si-nc embedded into a background with 
dielectric constant eo, an effective dielectric constant is 
defined as e e g = Eq + 4na/fl (fl is the structure vol- 
ume). £ e ff is a dielectric constant with SPEs included. In 
the case of noninteracting spheres, embedded into a back- 
ground with dielectric constant £o, the dielectric constant 
is given by^^ 

4esi - £o s 

£eff = £0 —77— 1 I 1 ) 

_£Si + 2e . 

where esi is the bulk Si dielectric constant. We refer to 
this way of calculating SPEs (that leads to Eq. ([1]) for 
spheres) as to the classical model (CM), in that QCEs are 
neglected, and classical electrostatics is used. The total 
dielectric constant e m of the Si-nc plus the background 
mixture is obtained by a weighted average between Eq. 
([T]) and the background dielectric constant. Introducing 
the Si-nc filling ratio /, the average is simply given by 

Em = /Eeff+(l-/)eo. (2) 

With the explicit substitution of e e ff, Eq. {2J gives the 
linearized Maxwell-Garnett effective medium equation. 12 
An enhancement of the model is obtained when a size de- 
pendent dielectric constant is used as esi- We call this 
latter the semiclassical model (SCM), in that, even if 
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SPEs are calculated classically, the model takes into ac- 
count QCEs. We stress here that esi does not include 
SPEs, and cannot be directly compared to experiments, 
unless SPEs are negligible. 

The available quantity in experiments is e m , the mix- 
ture dielectric constant. Therefore the comparison with 
the experimental data can be done following three steps. 
(1) The starting point is a dielectric constant, which ei- 
ther includes (SCM) or not (CM) the QCEs. (2) SPEs 
are introduced in a classical way, using Eq. (fTJ) in the 
case of spherical non interacting Si-nc. (3) A weighted 
average with the background dielectric constant is done, 
introducing the Si-nc filling factor. Points (2) and (3) 
constitute the basis of effective medium models, that in 
the case of noninteracting particles (low filling factors), 
reduce to the Maxwell- Garnett equation. Equation ([TJ) 
can be generalized to a frequency dependent dielectric 
constant. In this case esi is a complex quantity, and 
from e c ff the absorption cross section can be derived. 
The whole procedure goes under the name of dielectric 
model, since the materials are treated as continuous di- 
electric media. The model is commonly used to compare 
theoretical results and experimental dataJ^ The model 
is expected to fail either for small or for nanocrystals 
with a complex shape. For this reason it is important 
to calculate SPEs within an atomistic quantum mechan- 
ical framework, that reduces to the dielectric model for 
large size, regular shaped Si-nc. Moreover, an atomistic 
scheme is absolutely needed for the study of the effects on 
the optical response due to both the presence of inhomo- 
geneities on the surface, and the chemical nature of the 
passivantSfSili In the following we describe how SPEs are 
calculated using quantum mechanical atomistic theories. 

For years, dielectric constant calculations have been 
performed within an independent particle scheme, using 
the so called random phase approximation (RPA).— £ In- 
dependent particle RPA takes into account the QCEs, 
but it neglects the electron-hole interaction-^ and espe- 
cially LFEs. These are related to the inhomogeneous 
nature of the electron density induced by an external 
field. While in crystalline Si LFEs are basically due to 
the polarization of Si-Si bonds giving a negligible contri- 
bution to the dielectric constant, in Si-nc they are very 
important due to the presence of an interface between 
each nanocrystal and the background^ In such a case, 
LFEs are mostly due to SPEs, as it has been recently 
shown in the case of silicon quantum wires*! 7 - Due to the 
huge computational task, only recently LFEs have been 
calculated using first principles schemes i 14 i 17 ' 18 

We use a semiempirical tight binding scheme which 
gives the unique opportunity of studying structures made 
of hundreds of Si atoms. Details on the method and 
the approximation used for the position matrix elements 
can be found elsewherei 7 -^ It is just worth mention- 
ing that we make use of an sp 3 third nearest-neighbor 
parametrization, which provides for bulk silicon a good 
description of both the band structure and of the optical 
response. 7 Using the diagonal approximation for the posi- 
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FIG. 1: (Color online) Static dielectric constant calculated for 
a set of increasing size Si-nc using different levels of approx- 
imation: RPA, blue crosses; RPA+LF, black squares; semi- 
classical model, red plus symbols; classical model, green solid 
line. The lines are guides for the eyes. 

tion matrix elements (see Ref. and references therein) , 
we obtain the following expression for the frequency de- 
pendent dielectric tensor^ 

where f2 is the Si-nc volume, Rf 's are the positions of 
atoms composing the structure, and 7 are Cartesian 
components, SV, is the real space screened polarization 
matrix. S is defined by S — Pe^ 1 (matrix multiplica- 
tion), where Pij and are the polarization and the di- 
electric matrix in a tight binding representation, that we 
calculate following Ref. 19. When S is substituted by P 
(that is, when the polarization is unscreened), Eq. ([3]) 
gives the standard independent particle RPA expression. 

We thus study the Si-nc optical response using one of 
the following approximations: (i) RPA (QCEs are taken 
into account, SPEs neglected), (ii) RPA+LF (QCEs and 
LFEs are calculated within the tight binding framework), 
(iii) SCM (QCEs are calculated using the RPA approxi- 
mation, SPEs using the dielectric model), (iv) CM (QCEs 
are neglected, esi is the bulk Si dielectric function, SPEs 
are calculated using the dielectric model). 

In Fig. [1] we show the static dielectric constant for 
spherical hydrogenated Si-nc, as a function of the size. 
RPA, RPA+LF, SCM, and CM results are compared. 
First we note that RPA+LF calculations (black squares) 
well agree with the SCM approach (red plus symbols), 
in which SPEs are included into the RPA dielectric con- 
stant using Eq. |T]). It is remarkable that the semiclas- 
sical model gives an extremely accurate result not only 
in the large size region, where it is expected to work 
well, but also for very small Si-nc. Second, when LFEs 
are taken into account, the dependence on size of the ef- 
fective dielectric constant becomes very weak. Figure Q] 
shows that for the RPA+LF calculation there is a rela- 
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FIG. 2: (Color online) Imaginary part of the dielectric con- 
stant for a set of increasing size Si- nc. RPA, blue cross sym- 
bols; RPA+LF, black thick line; SCM, red thin line. For 
large nanocrystals, quantum effects become negligible. CM 
results are represented in the bottom part of the right-hand 
panel. They are obtained by using the bulk Si RPA dielec- 
tric constant (red-thin line), or experimental data (green-plus 
symbols) . 



tive change of about 3% when the size changes from 1.2 
to 2.4 nm, against 20% of the RPA curve. Moreover, 
the convergence to the large size limit is much faster, 
since it is reached already for small Si-nc. This means 
that in experimental situations in which SPEs are im- 
portant, SPEs play a role which is the opposite of QCEs, 
in that their contribution is size independent. For large 
size structures (dimensions greater than the Si exciton 
Bohr radius), QCEs are negligible and both RPA+LF 
and SCM tend to the CM limit. In Fig. [TJ for the CM 
calculation, the experimental value esi = 11.4 is used in 
Eq. fl}. A remarkable information emerging from Fig. 
Q] is the different behavior of RPA and RPA+LF results 
upon increasing the nanocrystal size. While RPA tends 
to the bulk Si dielectric constant, RPA+LF does not. 
This means that LFEs (that is SPEs) lead to a signifi- 
cant suppression of the static dielectric constant for each 
nanocrystal size. This is due to the fact that, at variance 
with the microscopic QCEs, SPEs give a macroscopic 
contribution. For this reason we could say that the effec- 
tive dielectric constant is a physically different property 
than the usual (RPA) Si-nc dielectric constant, and it 
is necessary to distinguish between calculations in which 
SPEs are either neglected or included. It is remarkable 
that SPEs comes out directly from the standard linear 
response theory, and there is no need of any further step. 

Figure [2] shows the imaginary part of the frequency 
dependent effective dielectric constant for several Si-nc 
and the classical model. A comparison between RPA, 
RPA+LF, and SCM is done on increasing the size. CM 
curves are calculated using the generalization to a disper- 
sive medium of Eq. ((T|) , in which esi is either the RPA 
(rod thin line) or the experimental data (green plus sym- 
bols) . The agreement between the two curves is good. As 
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FIG. 3: (Color online) Extinction coefficient of a d=1.8 nm 
Si-nc. RPA (blue dashed line) and RPA+LF (black solid line) 
are compared with the experimental data, taken from Ref. 2i 
(green symbols). 



discussed above for the static dielectric constant, LFEs 
are mostly due to SPEs, and the agreement between 
RPA+LF and SCM improves on increasing the Si-nc size. 
Furthermore, the convergence with the size is quite fast, 
and the Sii47Hioo, SiigiHug, and bulk Si spectra have 
the same trend in the low energy range, the only rele- 
vant difference being an overall rcdshift due to QCEs. 
From Fig. [2] it comes out that SPEs lead to a strong sup- 
pression of absorption in the visible range. It is worth 
stressing that the Si-nc optical absorption dramatically 
change when SPEs are taken into account. While RPA 
curves have structures quite similar to the bulk Si ab- 
sorption spectrum, with several peaks due to the most 
relevant interband transition energies, 7 RPA+LF curves 
have a monotonically increasing behavior in the whole 
energy range considered here. 

This behavior has been nicely confirmed by experi- 
ments. Figure [3] shows a comparison between the mea- 
sured and the computed extinction coefficient of a Si nc 
with a size of 1.8 nm (green symbols). We have already 
compared the RPA calculation (blue dashed curve) with 
this experimental datum! 7 - Here we discuss the role of lo- 
cal fields, by a direct comparison of the previous results 
with the RPA+LF curve, shown by a black solid line in 
Fig. [3J From the comparison, we can see that while the 
RPA curve decreases to zero after a main peak at about 5 
eV, the experimental curve shows an increasing behavior 
with a good agreement with the RPA+LF curve. 

It is important to remark that SPEs play a key role for 
the physical understanding of the spectra calculated by 
the time dependent density functional theory (TDDFT). 
When TDDFT absorption spectra were calculated for Si- 
nc, a strong discrepancy with respect to the DFT result 
came out^ Recently, an interpretation of TDDFT re- 
sults based on the classical theory has been giveni 23 i 24 
The DFT absorption spectra have peaks centered on the 
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FIG. 4: (Color online) Absorption cross section of Si3sH36. 
RPA+LF results, black thick line; semiclassical, red thin line; 
TDLDA results, green symbols. TDLDA results are taken 
from Ref. ED. 



energies relative to the most relevant filled-to-empty state 
electronic transitions, and they are in nice agreement 
with the RPA calculations performed with semiempirical 
methods. 2 -- The situation is very different when TDDFT 
calculations are performed. Si-nc acquire a featureless 



spectrum in the low energy range (0-6 eV), showing a 
nearly monotonically increasing behavior. 22 Recently, for 
small Si-nc it has been shown that TDDFT gives results 
similar to RPA+LF . 14 i 18 In Fig 0] a comparison of the 
present RPA+LF results with the TDLDA calculations 
of Ref. l2l| is shown. By considering the huge difference 
between the two schemes, the agreement is fair. 

In this paper, local field effects in the optical response 
of silicon nanocrystals have been studied. They have 
been shown to dramatically influence both the absorp- 
tion spectra and the static dielectric constant. The cause 
lies in the surface polarization effects that is a key ingre- 
dient in the study of finite systems i 2 ^ We have shown 
that the dielectric model well reproduces the atomistic 
results. Moreover, we have shown that the comparison 
with experiments and first principles results is good. 
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